modus <- function(v) {
  
  anzahl = max(v)
  m <- matrix(0, anzahl, 1)
  
  for (i in (1:length(v))) m[v[i],1] <- m[v[i],1] + 1
  
  ausgabe <- c()
  
  for (i in (1:length(m))) {
    if (m[i,1] == max(m)) ausgabe <- append(ausgabe,i)
  }
  
  return (ausgabe)
}


v <- c(6,8,14,21,25,14,5,6,2,1,3,6,16,20,28,18,4)

print("Analyse der Kaffeedaten")
print(paste("Mittelwert:", mean(v)) )
print(paste("Median:", median(v)) )
print(paste("Modus: ",modus(v)))

